home *** CD-ROM | disk | FTP | other *** search
/ Collection of Tools & Utilities / Collection of Tools and Utilities.iso / fortran / f2c_src.zip / F2C / LIBF77 / POW_ZZ.C < prev    next >
Encoding:
C/C++ Source or Header  |  1991-06-10  |  329 b   |  18 lines
  1. #include "f2c.h"
  2.  
  3. VOID pow_zz(r,a,b)
  4. doublecomplex *r, *a, *b;
  5. {
  6. double logr, logi, x, y;
  7. double log(), exp(), cos(), sin(), atan2(), cabs();
  8.  
  9. logr = log( cabs(a->r, a->i) );
  10. logi = atan2(a->i, a->r);
  11.  
  12. x = exp( logr * b->r - logi * b->i );
  13. y = logr * b->i + logi * b->r;
  14.  
  15. r->r = x * cos(y);
  16. r->i = x * sin(y);
  17. }
  18.